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C^) ' Abstract. We revisit the issue of ripples in the primordial power spectra caused by 

0O , trans-Planckian physics, and the potential for their detection by future cosmological 

probes. We find that for reasonably large values of the first slow-roll parameter e 

( > 0.001), a positive detection of trans-Planckian ripples can be made even if the 

K^ , amplitude is as low as 10~ 4 . Data from the Large Synoptic Survey Telescope (LSST) 

^ ' and the proposed future 21 cm survey with the Fast Fourier Transform Telescope 

(FFTT) will be particularly useful in this regard. If the scale of inflation is close to 
its present upper bound, a scale of new physics as high as ~ 0.2 Mp could lead to 
observable signatures. 
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1. Introduction 

Observations of the temperature and polarisation anisotropies in the cosmic microwave 
background (CMB) radiation strongly support the idea that the initial conditions for 
structure formation were set by an earlier period of cosmic inflation [1-4]. During this 
period, the rapid expansion of space caused quantum fluctuations of the inflaton field to 
be blown up to macroscopic scales. The initial state of these quantum fluctuations 
is typically assumed to be given by the Bunch-Davies vacuum of de Sitter space 
at a time when k ^> aH, i.e., the physical wavenumber k/a is much larger than 
the Hubble expansion rate. Along with the usual assumptions of single-field, slow- 
roll inflation, this leads to the almost scale-invariant spectrum of Gaussian adiabatic 
scalar perturbations that constitutes the "primordial" state of perturbations in today's 
cosmological concordance model [5]. 

However, in an expanding space-time the notion of a vacuum is not unique. This can 
easily be understood by observing that the expanding space-time makes the Hamiltonian 
depend explicitly on time, so there is no time-independent lowest energy state that can 
serve as a vacuum at all times. We might then worry about how sensitive our predictions 
are to the choice of initial state. The choice of vacua other than the Bunch-Davies 
one leads in principle to unacceptable ultraviolet divergences in the observables [6], 
although this deficiency can be cured by introducing a high-energy cutoff to the theory. 
In practice, such a cutoff will be given by the Planck scale, or possibly even lower, by 
a scale of new physics A (if the cutoff is too close to the Planck scale, an unorthodox 
vacuum choice may in fact lead to backreaction problems [7-12]). An ultraviolet cutoff 
in the inflaton fluctuation modes violates energy conservation and has to be associated 
with a source term for the fluctuations, in order to account for the modes redshifting 
across the new physics hyper surf ace [10]. 

Unless inflation started only just before the wavelengths corresponding to today's 
observable scales left the horizon, these fluctuations can be mapped to wavenumbers 
larger than A during inflation and have emerged from above the new physics hyper- 
surface. It is thus very well conceivable that the fluctuations bear an imprint of 
the unknown new physics [13,14]. In particular, the modes might be created in a 
non-standard vacuum state at the new physics hypersurface. Danielsson proposed a 
prescription of how to construct this initial state by minimising the Hamiltonian on 
this hypersurface [15]. However, it was found in references [16,17] that the state 
constructed in this fashion is somewhat ambiguous, since it is not invariant under 
canonical transformations of the Hamiltonian. In fact, it has been argued from a purely 
effective field theory point of view that the effect on the primordial spectrum should be 
smaller than that proposed by Danielsson [18]. Eventually, only a better understanding 
of the thinning of ultraviolet degrees of freedom in quantum gravity [19] can help us 
solve the puzzle of the correct choice of initial conditions and how the fluctuations 
emerge during inflation, as in the proposed trans- Planckian damping mechanism [20] 
(see also [9]). 
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Nonetheless, the generic signature of choosing a non-Bunch-Davies vacuum state 
appears to be a superimposed oscillation of the primordial spectra ("ripples"), whose 
amplitude is suppressed by a power of H/A. For definiteness we will consider 
Danielsson's model as a case study in the following. 

The question of whether traces of these oscillations can possibly be detected in 
present or future CMB anisotropy data has been the subject of a number of recent 
studies [21-30]. While few will disagree that the CMB anisotropics are at present the 
most powerful probe of the primordial power spectrum, the advent of high-redshift 
surveys of the large scale structure (LSS) distribution in the not- so- distant future may 
help to enlarge the detection window. Planned galaxy redshift surveys such as the 
Large Synoptic Survey Telescope (LSST) [31], or surveys of the distribution of neutral 
hydrogen using the 21 cm spin-flip line (e.g., the Murchison Widefield Array (MWA) 
[32] , the Square Kilometre Array [33] , or the Fast Fourier Transform Telescope (FFTT) 
[34]) have the advantage of mapping out the density perturbations in three dimensions 
over large volumes. They are less affected by sampling variance (or cosmic variance), 
which is ultimately the limiting factor for CMB data over a large range of scales. In this 
work, we investigate how a combination of these probes with CMB data can enhance our 
potential to discover the footprints of new physics beyond the energy scale of inflation. 

We shall begin by refreshing the reader's memory of the Danielsson model in 
section 2. Some general considerations about the phenomenology of oscillating spectra 
will be presented in section 3, before we turn to the questions of under which conditions 
their traces can be detected in CMB data (section 4), or, if no detection is possible, 
what constraints can be placed on the parameters of the model (section 5). In sections 6 
and 7 we discuss how future galaxy redshift surveys and 21 cm surveys can contribute 
to improving the chances for a detection of ripples. We conclude in section 8. 

2. The model 

The spectrum of CMB perturbations predicted from inflation is usually expressed in 
terms of the spectrum of co-moving curvature perturbations, since this is a conserved 
quantity on super-Hubble scales in single field slow-roll inflation. In the spatially flat 
gauge, the curvature perturbation TZ(k) can be related to the fluctuation of the inflaton 
field 5(f>(k) by the relation S(j)(k) = (<fi/H)7Z(k), where H is the Hubble scale and is 
the classical expectation value of the inflaton field. 

Let us remind the reader of the Danielsson initial condition and at the same time 
generalise the derivation of the spectra to the case of standard slow-roll inflation. 
Since there has been some confusion about the correct generalisation to slow-roll in 
the literature, we repeat some of the details. 

In the spatially flat gauge we can write the metric in the form 

ds 2 = a 2 (T)(dT 2 -dx 2 ), (2.1) 

and after transforming to the canonical variable /i = a5<ft, the perturbation of the inflaton 
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field satisfies the Mukhanov equation, which takes the form of the equation of motion 
of a harmonic oscillator with a time- dependent mass [5] , 
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where v = 3/2 + 3e — 77, and e and rj are the usual slow-roll parameters which we assume 
to be constant. The conjugate momentum to fi k is given by [5] 
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where z is defined as z = v 2ea. 

Like for the ordinary harmonic oscillator, it is convenient to quantise the system by 
promoting the fields to operators and expanding them in terms of raising and lowering 
operators 
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where the vacuum is defined as the state annihilated by the lowering operator 

h (r) |0) T = 0. (2.5) 

Now, the vacuum at any later time is determined by a Bogoliubov transformation 

a k (r) = a k (r) a k (r ) + f3 k (r) aL k (r ) , 



(r) = a^ (r) a k (r ) + ^ (r) a_ k (r ) 



(2.6) 



where a^ and /3fc are the Bogoliubov coefficients satisfying the normalisation condition 
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in order for the commutation relations to be preserved in time. The solution to the 
dynamical equations is most easily found by defining a set of mode functions 
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h (V) = ^r K (r) + # (r)] , 
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which satisfies the classical equation of motion, as can be seen by employing the 
Heisenberg equation of motion. One can then invert these relations to find 

to 0?) = /* (v) Ofe {m) + fk (v) a~ k (m) > 

TTfc (^) = 9k (v) ^k (Vo) + 9* k (v) a -k (Vo) ■ (2.9) 

By minimising the Hamiltonian written in terms of the variable <p and its conjugate 
momenta (note that minimising the Hamiltonian written in terms of the Mukhanov 
variable gives a different result [16,17]), one can show that the instantaneous vacuum 
|0) ro that minimises the energy at r = r is determined by 

9k (r ) = -ikf k (r ) . (2.10) 
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This is equivalent to the zeroth order adiabatic vacuum, and the condition is equivalent 
to requiring /3(r ) = 0, which one can use to normalise the general solution to the 
Mukhanov equation 
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where A k and B k are constants of integration to be fixed by the physical initial conditions 
and £ = —{kv + \)n. Solving the equations for the Bogoliubov coefficients a k and (3k, 
and expanding in 1/x, one obtains 
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As mentioned above, the Danielsson choice of vacuum requires that we put (3 k {tq) = 
at some initial moment Tq. This implies that 
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One can then calculate the spectrum of curvature perturbations P-r = (l/2e)P ( / ) , 

which becomes 
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where we have taken the super-horizon limit t — > and evaluated the expression at 
horizon crossing k = aH , using r = — l/[aH (1 — e)] in the prefactor. Inserting \ gives, 
up to a phase, 
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where x = —krj . 

This expression agrees with the straightforward generalisation of the corresponding 
calculation for power-law inflation in references [16, 17], but disagrees slightly with the 
result of [35], who found that the expression does not depend on rj. This can be traced 
to an inaccurate approximation.^: 

Following the proposal of Danielsson, we should fix the initial condition of a given 
mode when that mode crosses a fixed cutoff, i.e., 

\x \ = A/ H , (2.19) 

where factors of order one can be absorbed in the arbitrary cutoff A. Here H is the 
Hubble parameter at r = r , when the mode crosses the new physics hypersurface 
defined by the ultraviolet cutoff A. Since H is not the same for all modes because of 
the slowly decreasing H, this will introduce /c-dependent oscillations in the spectrum. 
It can be shown that this leads to xq oc k e , since Hq oc k~ e [35]. One can then normalise 
H to H to obtain 
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where £ = H /A. The phase 
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(2.21) 



needs to be added as a free parameter since any 
dependence of the final result on the (arbitrary) choice of pivot scale fco is non-physical 
and must necessarily be spurious. A similar expression exists for the tensor power 
spectrum, i.e., P h (k) = Pj? D {k)(- ■ ■). 



3. Preliminary considerations 

In the limit eCl, the argument of the sine function in equation (2.21) can be expanded 
in powers of e. To the lowest order in e and rj, 



P**(k) ~ P£Wk) \l + t 
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sin 




(3.1) 



i.e., the oscillations are in lnfc, with frequency uj = 2e/£, or equivalently, oscillation 
length L = ir£,/e, and we have defined 0' = <fi + 2(1 + e)/£. 

| The procedure outlined in the footnote on page 7 of [35] is slightly misleading. Although the zeroth 
order adiabatic approximation is not exact, its solution is only used as a boundary condition to normalise 
the exact solution at the initial time giving the zeroth order adiabatic vacuum solution, which is an 
exact solution. In the end we want to end up with an exact solution, otherwise the effect would just 
be an effect of an inaccurate approximation. 
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3.1. Constraints on the oscillation length 

The totality of CMB and LSS observations spans roughly four orders of magnitude in k- 
space between /c min ~ 1CT 5 h Mpc -1 and fc max ~ 1CT 1 h Mpc -1 . In natural log space this 
range is approximately hi(/c max //c min ) ~ ln(10 4 ) ~ 9.2. For the individual experiment 
we take the observable k range to be that in which the measurements are effectively 
noise-free and limited only by cosmic variance. 

For a CMB survey that is cosmic variance-limited up to £ = 2000 (the "CVL 
survey"), we have /c max — 0.15 h Mpc -1 within the concordance ACDM framework. The 
Planck satellite also has similar specifications in its TT measurements [36]. The WMAP 
TT power spectrum is cosmic variance-limited only up to I ~ 500 [37], corresponding 
to /c max ~ 5 x 10 -2 h Mpc -1 . Therefore, 

f 6 9 WMAP 
i?CMB = ln( W*W) - | g Q ; pknck ^ CVL5 (3.2) 

where in both cases we have used k min ~5x 10 -5 h Mpc -1 . 

If trans-Planckian corrections to the primordial power spectrum are to be observed 
as an oscillatory signal, we must be able to fit in at least three nodes inside the observable 
range -Rcmb- Anything less will not qualify as an oscillation, and the corrections will 
at most show up as a spectral tilt or perhaps a running, or in the worst case, as a 
normalisation. Depending on the phase <fi three nodes is equivalent to (1 -j-5/4)L, where 
L is the oscillation length defined immediately after equation (3.1). This then sets an 
upper limit on the range of £/e, 

£ < 4 f 1.8, WMAP, 

7 <-i?CMB-j 205 Planck and CVL) ( 3 - 3 ) 

in which some vaguely oscillatory effect may be observable. 

A lower bound on £/e arises from the finite widths of the CMB "window functions" . 
The CMB anisotropy spectra can be written schematically as 

C 1 / = £ fd\nkT e ' a (k)Tt(k)P a (k), (3.4) 

where i,j = T,E,B, and a sums over scalar and tensor modes. If the width 
of T\ ,0i (k)Tp a (k) is much larger than the oscillation length L, then we expect the 
oscillations to be wiped out, or, at least, suppressed by a factor of a few. For the 
ACDM model, the log width oiT l e ' a (k)T^ a {k) at £ ~ 1000 is ~ 0.07, and increases with 
decreasing £. Therefore, in order to discern an oscillatory signature, we must have 

-> 0.022. (3.5) 

We note that the window functions for the LSS power spectrum from present galaxy 
redshift surveys are even broader. For the power spectrum of the Sloan Digital Sky 
Survey Luminous Red Galaxy (SDSS LRG) sample published in reference [38], the 
window functions have a log width of ~ 0.33, leading to £/e > 0.1. This means that 
present galaxy redshift surveys are not very useful for resolving small ripples. 
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3.2. Constraints on the oscillation amplitude 

Suppose the true primordial power spectrum does have trans-Planckian ripples. One 
can then ask, how large does £ need to be so that a smooth spectrum is a "bad" fit to 
the data. This can be estimated as follows. 

The effective y 2 for a cosmic variance-limited CMB experiment in one mode (e.g., 
TT) is given by 

X e 2 ff = / sky E(2^+l) R^ + ln^ > ( 3 - 6 ) 

1=2 V G ^ U t / 

where C£ bs denotes the observed CMB anisotropy spectrum, Cg the theoretical 
predictions, and we have included a factor / s k y — 0.7 to account for the inevitable 
sky cut. The difference in xls between the smooth and the oscillatory models is then 

Ch / /^<obs /^obs /^smooth \ 

Ax^ ff = /sky E (2^+1) 7^-7^ + 1*-^- . (3-7) 

where C| mooth and C£ sc denote the theoretical predictions of the two respective models. 
We express the oscillatory anisotropy spectrum C^ sc in terms of the smooth one as 

Cf c = Cf ooath (l + A t ), (3.8) 

where A^ <C 1 contains the oscillatory features and in general must be calculated by 
folding in the appropriate transfer functions with the primordial power spectrum (2.21). 
Here, we estimate its effects by substituting A rms for A^, where 

is the rms amplitude. Then, keeping only terms to the lowest order in £, we get 

£2 cmax 

AXeff - -T /sky E U + 1 
4 1=2 

= "j" /sky (/max + 3J(c max — 1) — ~Tjskyt max , (3.10) 



where we have also used the expectation value (C^ bs ) = C^ sc . 

As said, the expression (3.10) assumes we are measuring in one mode only (e.g., 
TT). To account, approximately, for additional modes (e.g., TE, EE), we need only 
to sum their contributions to A%^ ff . We note that this procedure is strictly not correct 
since the TT, TE and EE modes in principle correlated, and a simple summation would 
tend to overestimate A%g ff . However, since we are interested only in order magnitude 
estimates, a small overestimation is not a major concern. 

Suppose our TT measurements are cosmic variance-limited up to £ max = 2000 (e.g, 
Planck). Then a "2<r" detection of trans-Planckian ripples, i.e., Ax^ ff > 4, would require 

£> 0.0024. (3.11) 

If in addition TE and EE are cosmic variance-limited up to the same ^ max , 

£> 0.0014. (3.12) 
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Figure 1. Schematic showing the observability of trans-Planckian ripples. The 
blue/solid box on the right represents the region in which the oscillation length is so 
long that £ becomes degenerate with other parameters describing the primordial power 
spectrum. The blue/solid box on the left denotes the region in which the ripples are 
wiped out by the window functions of the experiment because of too short an oscillation 
length. The red/hatched box on the left encompasses the region in which the oscillation 
amplitude is too small compared with the sampling error (i.e., cosmic variance). Two 
vertical dashed lines demarcate the only window in which trans-Plankian ripples may 
be detectable. 



The corresponding limit for a WMAP-like measurement is 
£ > 0-0096, 



(3.13) 



noting that the WMAP TT measurement is effectively noise-free only up to £ max ~ 500. 
Finally, in deriving these amplitude limits we have assumed that £ is not degenerate 
with other cosmological parameters. This is a good approximation if £ satisfies the 
bounds (3.3) and (3.5) on the oscillation length. In other words, if trans-Planckian 
effects manifest themselves as oscillations in the observable range, then no other 
parameter describing the primordial power spectrum can mimic their effects on the 
cosmological observables. Trans-Planckian ripples might conceivably be confused with 
baryon acoustic oscillations (see reference [39] for a related scenario where this is the 
case). We do not, however, find this to be a problem since the ripples are oscillations 
in ln/c, while the baryon acoustic oscillations are oscillations in k. Results from our 
numerical parameter error forecasts corroborate this claim. 



4. Can we detect ripples? 

A simple picture emerges from the above considerations (figure 1). For a given e, 
oscillation length requirements define a window in £ in which trans-Planckian ripples are 
potentially observable. Amplitude requirements define a second, e-independent window. 
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Table 1. Ranges of £ in which trans-Planckian effects can be observed as oscillations 
for WMAP, Planck, and CVL for fixed values of e. The parentheses denote those cases 
that are in reality unobservable, either because of an undetectably small e, or because 
the expected measurement error on £ exceeds the width of the detection window. 

e WMAP Planck CVL 

0.02 (0.0096 < £ < 0.036) 0.0024 < £ < 0.041 0.0014 < £ < 0.041 

0.01 (0.0096 <£<0.018) 0.0024 < £ < 0.020 0.0014 < £ < 0.020 
0.001 (0.0014 <£ < 0.0020) 

0.0001 



Only when the two windows overlap can one expect to see ripples. Table 1 summarises 
the £ ranges in which ripples can be observed for WMAP, Planck, and a CMB survey 
that is cosmic variance-limited in TT, TE and EE up to £ max = 2000 (CVL) for various 
fixed values of e. 

Clearly, a crucial prerequisite for ripples detection is that e must be sizeable. With 
e fixed at a value as large as > 0.01, WMAP, Planck, and CVL all have the potential 
to see trans-Planckian ripples, although the observable window for WMAP is so small 
that taking into account the expected measurement error on £, 

A£~£ min , (4.1) 

where £ m i n denotes the minimum detectable £, would most likely close it. If we fix 
e = 0.001, ripples are unobservable for WMAP and Planck, but are in principle still in 
range for CVL (barring the finite A£). No detection window exists for e = 0.0001 or 
smaller, since the associated oscillation lengths are so large that £ always falls in the 
region degenerate with the spectral tilt and/or the normalisation. 

The question then is, how well do we expect to measure e from CMB data, since 
by the same arguments any measurement of e that is consistent with e = means that 
no window exists for the detection of trans-Planckian ripples. A null detection of e 
by WMAP, e < 0.02 (95% C.L.), essentially rules out the observation of ripples by the 
current generation of CMB probes. The projected la sensitivities to e for Planck and 
CVL from our parameter error forecasts (see Appendix A for details) are Ae ~ 0.005 
and Ae ~ 0.0015 respectively for an undetectable £. Some variations exist if there is a 
detection of £ (empirically, Ae/e ~ A£/£). In table 1, we enclose in parentheses those 
scenarios that are most likely unobservable either because of an undetectably small e, 
or because the expected error on £ exceeds the width of the detection window. 

To confirm our analytic estimates we perform several rigorous parameter error 
forecasts. Figure 2 shows the A%^ ff curves for three different fiducial models, assuming 
the CVL experiment. Technical details of the forecast and definitions of the statistical 
quantities can be found in Appendix A. 

For fiducial model (i) with ^ = 0.01 and €m = 0.01, the prospects for ripples 
detection are excellent. An e fid so large can be pinned down by CVL at better than 5cr. 
This enables a detection of £ with a correspondingly high statistical significance: between 
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Figure 2. Axis curves as functions of £ and e for three different fiducial models for 
the CVL experiment: (i) £ fid = 0.01, e fid = 0.01 (red/solid), (ii) £ fid = 0.01, e fid = 0.001 
(blue/dashed), and (iii) £g d = 0.001, efi d = 0.01 (green/dotted). The three horizontal 
lines indicate, from top to bottom, A^g = 9,4, 1 (i.e., 3a, 2a, la). 



the best-fit £ and at £ — > (i.e., a smooth spectrum), we find Axl s ~ 40. However, if e^ 
is lowered to 0.001 (model (ii)), no detection of £ is possible, since e = 0.001±0.0015 (la) 
is consistent with e = so that ^ = 0.01 lies in the degenerate region. 

Model (iii) with ^ = 0.001 and egd = 0.01 is also futile, even though the projected 
la error on e, Ae ~ 0.0015, makes e^ well within experimental reach. The reason here 
is that £fid = 0.001 leads to an unobservably small oscillation amplitude. Interestingly, 
although the £m of model (iii) is beyond CVL's reach, the very large and detectable 
egd still makes it possible to exclude a range of £ values from observations. This is in 
contrast with model (ii), for which an undetectably small e&& means that no conclusions 
can be drawn about £ (to be discussed in more detail in section 5). 

The results of figure 2 are in excellent agreement with our analytic estimates 
presented in table 1. Using the same arguments we can also conclude that if £m = 0.01 is 
to be detected by CVL, the underlying model must have e^ > 0.005. This result is fully 
consistent with a similar claim by Easther, Kinney and Peiris (EKP) [27], who found 
that £ EKP = 0.01 is detectable at 2a by a CMB experiment that is cosmic variance- 
limited in TT, TE, EE and BB up to £ max = 1500 only if r > 0.1 (or e > 0.006). § 

§ EKP used a power spectrum slightly different from our equation (2.21): P(k) = P BD (k)(- ■ -) 1 / 4 . 
This means our £/e constraints from oscillation length considerations apply also to £ EKP /e, while our 
lower limits on £ from amplitude arguments need to be multiplied by a factor of four for use with £ . 
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5. Can we constrain ripples? 

Suppose we detect e = 0.01, and the true primordial power spectrum has no trans- 
Planckian ripples. Then using the same amplitude arguments as those in section 3.2, 
we can put a la constraint on £ of 

i< 0.0014 (5.1) 

for the CVL experiment. If we take instead the EKP ripples and experimental set-up, 
the same procedure gives a la limit of £ EKP < 0.0064, which is in good agreement with 
EKP's £ EKP < 0.004 for a scenario with r = 0.15. 

However, this upper limit is not the whole story. This is because if £ is very large 
(£ ^> e), then for any given e<l, there exists an exact degeneracy between £ and 0. 
To see this we perform a small x = (2e/£) ln(/c//c ) expansion of the sine function in 
equation (3.1). This gives us a set of effective parameters 

Af t = A s (l + tsm<f>'), 

C ff , e(2cos0'-£sin0') 

n sA = n s,t + - 



af = a s 4 



1 + £ sin (f)' 
2e(n s — 1) (2 cos 0' — £ sin < 



S * ' IXfclr,^ 



1 + £ sin < 

eff , 2era i (2cos0'-£sin0') 

a t = a t + . • (5.2) 

1 + £ sm 0' 

Applying the consistency relations A^ /Af = A t /A s = 16e and rff = rit = — 2e, we see 
that for whatever £ and e we choose, setting 0' = arctan(2/£) always returns the exact 
solutions A s>t = Afl, n s> t = nf, and a S)t = off t . 



More (approximate) solutions exist if we ignore the consistency constraints on n 



off 

t 
and a* (a good approximation since future data are unlikely to resolve the tensor 

power spectrum) , or if the running of the scalar spectral tilt a s is treated also as a free 

parameter. Thus, unless we have a reason to fix (0, n s ,A s ) to some particular values, 

it is always possible to find a fit in the large £ S> e limit that is almost or even exactly 

as good as that of a smooth spectrum. For this reason a lower limit on large £ should 

exist in addition to the upper limit (5.1) on small £. 

Unfortunately there is no simple way to estimate this lower limit accurately, since 

before we enter the exactly degenerate region at £ ^> e, there exists a region at £ ~ e 

in which £ already exhibits some (partial) degeneracy with other parameters describing 

the primordial power spectrum. The small x expansion limits the exactly degenerate 

region to x < 1, or, 

yi?CMB<l, => £>18e (5.3) 

for the CVL experiment. A less stringent limit can be set by the same oscillation length 
arguments that led to equation (3.3), i.e., 

£ > 2.0c, (5.4) 
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Figure 3. Same as figure 2, but for the models (i) £fid = 0.0, e^d = 0.01 (red/solid), 
and (ii) £ fid = 0.0, e fid = 0.001 (blue/dashed). 

again for CVL. The true limit likely lies somewhere in between (5.3) and (5.4). 

For e > 0.01, these lower limits on large £ are purely academic since the trans- 
Planckian spectrum used here assumes ^ < 1. They are important, however, in the 
case of small e, since the lower limit (5.3) or (5.4) may run into the upper limit (5.1), 
thereby closing the window in which some finite £ may be a "bad fit" to the data. For 
example, no window exists for the case of a fixed e = 0.0001. In the same vein, if e = 
is consistent with data, then we can draw no conclusions about £. 

Figure 3 shows the Ax^ ff curves for two fiducial models, (i) ^ = 0.0,e fid = 0.01, 
and (ii) ^ = 0.0, e^ = 0.001, assuming the CVL experiment. The bounds on £ for 
model (i) are in excellent agreement with the predictions (5.1) and (5.4). Model (ii) has 
an undetectably small ega- Correspondingly, no bounds can be set on £ in this case. 

Lastly, we note that EKP found an upper limit on £ even in the case of an 
undetectably small e, in apparent contradiction with our results. We believe this to be 
a consequence of the Bayesian statistics used in their analysis. In the small amplitude 
region at small £, all values of the phase <fi provide a good fit to the data, since the 
oscillations are in any case unobservable. In the degenerate region at larger £ values, 
however, only certain choices of give good fits. This reduces the volume of the posterior 
probability distribution in the direction at large £ values, which are subsequently 
disfavoured in a Bayesian analysis. 



6. Galaxy redshift surveys 

Galaxy redshift surveys offer at present the most reliable probe of the large scale 
structure distribution in the low-redshift universe (z < 1). Clustering statistics of 
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galaxies from the catalogues of the Two-Degree Field Galaxy Redshift Survey (2dFGRS) 
and the Sloan Digital Sky Survey (SDSS) are now routinely used to constrain 
cosmological parameters [38,40]. Many planned/proposed surveys will probe an even 
larger volume of the universe and at higher redshifts in the future. Here we consider 
the potential of the Large Synoptic Survey Telescope (LSST; first light 2015) [31]. 

The LSST is a full-sky survey looking out to a maximum redshift of z max ~ 3. 
Within the ACDM framework this corresponds to a survey volume of V ~ 100 h~ 3 Gpc 3 . 
Formally the noise-free region is defined as the range of k values for which the power 
spectrum signal P^ dominates over the Poisson shot-noise. We take this region to be 
k < k max ~ 0.2 h Mpc" 1 . We caution however that nonlinear evolution such as mode- 
coupling may yet suppress the trans-Planckian ripples or cause them to completely 
disappear beyond some k value smaller than k max . 

The largest scale accessible to the LSST corresponds approximately to k m \ n ~ 
27T/V 1 / 3 ~ 1.4 x 10 -3 h Mpc" 1 . The total k range probed by the survey therefore spans 
barely two orders of magnitude. To maximise the spectral coverage we combine the 
measurements from the LSST with those from a CMB experiment. Oscillation length 
arguments then impose an upper limit on £/e for which trans-Planckian ripples can be 
observed by CMB+LSST, 

^<^m(S T ASr)-2.1, (6.1) 

a number similar to the constraints (3.3) from the CVL experiment alone. 

The effective volume of the LSST is some 200 times larger than that of the SDSS 
LRG sample. Thus we can reasonably expect a corresponding factor of 200 1 / 3 ~ 6 
decrease in the widths of the fc-space window functions. The window functions of the 
SDSS LRG power spectrum has a log width of ~ 0.33, from which we derive a log 
width of ~ 0.055 for the LSST. This sets a lower limit on the oscillation length of the 
trans-Planckian ripples, or, equivalently, 

^> 0.018, (6.2) 

a limit comparable to the natural window functions of the CMB (3.5). 

To estimate the lowest observable £ from amplitude arguments we adopt the 
following expression for the effective \ 2 °f a sampling variance-limited LSS survey [41]: 



V f k -^JP k -Pu h& 



2 



If the true primordial power spectrum has trans-Planckian ripples, then the A%^ ff 
between a smooth and an oscillatory spectrum is 



> max 



AXeff = Jt~ 2 t~~ k 2 dkAl (6.4) 

where A& <^C 1 encapsulates the oscillatory part of the power spectrum. As in section 3.2, 
we substitute A*, with the rms amplitude A rms = £/v^2, thus leading to 

V £ 2 V £ 2 

AXcS ~ (27r)2lf ^ max ~ kmin ^ ~ j2Tiy~6 kniax ' ( ' 6 ' 5 ' ) 
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Table 2. Same as table 1, but for Planck+LSST, and CVL+LSST. 
e Planck+LSST CVL+LSST 



0.02 9.9 x 10~ 4 < £ < 0.042 8.5 x 10~ 4 < £ < 0.042 

0.01 9.9 x 10" 4 < £< 0.021 8.5 x 10~ 4 < £< 0.021 

0.001 (9.9 x 10~ 4 <£ < 0.0021) (8.5 x 10~ 4 < £ < 0.0021) 
0.0001 



Adding this Ax^ ff to its CMB counterpart (3.10) yields the total A%^ otal for the 
combination CMB+LSST. 

If we demand Axj? otal > 4 for a 2<r detection, then 

f 9.9 x 10" 4 , Planck + LSST, 
^~\8.5xl0- 4 , CVL + LSST. ^ ' ^ 

Compared with the corresponding limit for Planck alone (3.11), the projected constraint 
on £ for Planck+LSST is more than a factor of two better. Adding LSST to CVL also 
yields some improvement over the limit from CVL alone (3.12). Thus, with the LSST, we 
will be able to reach down to the £ ~ 0.001 level for the first time. Table 2 summarises 
the ranges of £ in which trans-Planckian ripples can be observed for Planck+LSST and 
CVL+LSST for various fixed values of e. 

The remaining question is how sensitive are these combinations of CMB and LSS 
probes to e. In the case of a non-detection of £, we find from our parameter error 
forecasts la errors of approximately Ae ~ 0.0015 and Ae ~ 0.005 for Planck+LSST 
and CVL+LSST respectively. These are essentially the same sensitivities reachable by 
Planck and CVL alone. In contrast, a detection of ^d = 0.01 by Planck+LSST or 
CVL+LSST does lead to some improvement in the sensitivity to e, Ae ~ 0.0007, for 
Cfid = 0.01. In either case, e = 0.001 remains out of reach. 

7. 21 cm surveys 

Recently there is a growing interest to map the three-dimensional distribution of 
neutral hydrogen via the observation of brightness temperature fluctuations in the 
21 cm spin-flip transition line. A number of radio arrays have been planned/proposed 
for this observation, targeting specifically the redshift range 6 < z < 20 [32-34]. At 
these redshifts spatial fluctuations of the 21 cm signal are in general determined by a 
combination of the underlying matter distribution and reionisation physics. However 
if the signal originates at a time when the neutral hydrogen spin temperature is much 
larger than the CMB temperature (i.e., z < 10) and if at the same time the intergalactic 
medium is completely neutral, the 21 cm spatial fluctuations are an exact tracer of the 
perturbations in the underlying matter density field up to a bias factor. It is within this 
very optimistic scenario that we consider the detectability of trans-Planckian ripples by 
21 cm surveys. See reference [42] for a review on the physics of the 21 cm transition. 
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Table 3. Same as table 1, but for Planck+MWA, Planck+SKA, and CVL+FFTT. 

e Planck+MWA Planck+SKA CVL+FFTT 

0.02 0.0024 < £ < 0.040 0.0014 <£< 0.048 7.9 x 10" 5 <£< 0.050 

0.01 0.0024 <£< 0.020 0.0014 <£< 0.024 7.9 x 10~ 5 <£< 0.025 

0.001 (0.0014 < £ < 0.0024) 7.9 x 10~ 5 < £ < 0.0025 

0.0001 (7.9 x 10~ 5 <£ < 2.5 x 10~ 4 ) 



Because they observe at high redshifts, 21 cm surveys tend to have large effective 
volumes and can in principle probe larger wavenumbers k of the LSS power spectrum 
without impediments from nonlinear evolution. However, these surveys are limited at 
the low k end of the spectrum because of residual foregrounds. Typically one expects the 
foreground to swamp out the 21 cm signal at k < 0.07 h Mpc -1 . Thus, it is necessary 
to combine 21 cm with CMB observations to optimise k m i n . 

To estimate /c max we use figure 3 of reference [43] which shows the effective volumes 
of the Murchison Widefield Array (MWA) [32], the Square Kilometre Array (SKA) [33], 
and the Fast Fourier Transform Telescope (FFTT) [34]. We define /c max as the point at 
which the effective volume drops below one half of the maximum effective volume for each 
survey (effectively where the thermal noise of the radio telescope becomes comparable 
to the 21 cm signal). For MWA, SKA and FFTT, we find & max ~ 0.15,0.6, 1 h Mpc" 1 
respectively. Oscillation length arguments then lead to the bounds 

f 2.0, CMB + MWA, 
i < — ln^aT ASD *\ 2.4, CMB + SKA, (7.1) 

[ 2.5, CMB + FFTT. 

To estimate the lowest observable £ by our 21 cm surveys we use again the Axis 
expression (6.5), but replace V with the effective volumes from figure 3 of reference [43] 
(i.e., V ~ 1.5 h~ 3 Gpc 3 for MWA and SKA, and V ~ 150 h~ 3 Gpc 3 for FFTT). 
Again, we compute Ax^ otal by adding Axlg to its CMB counterpart. Demanding that 

A Xt 2 otai Z 4 then g ives 

0.0024, Planck + MWA, 

0.0014, Planck + SKA, (7.2) 

7.9 x 10" 5 , CVL + FFTT. 

Table 3 contains a summary of the £ ranges detectable by Planck+MWA, Planck+SKA, 
and CVL+FFTT. 

With a potential to reach down to £ < 10~ 4 , CVL+FFTT improves the sensitivity 
of CVL alone by a factor of 20. That CVL and FFTT observe at different wavenumbers 
(but with some overlap) also allows their combination to probe, for some given e, larger 
values of £ than either experiment alone. Of particular note is that the improved limits 
will create, for the first time, a sizeable detection window for e values as small as 0.001. 
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8. Conclusions 

We have investigated the sensitivity of various future cosmological probes to "ripples" 
in the primordial power spectrum due to "trans-Planckian" physics. The magnitude of 
this effect is determined by £ = Ho/ A, i.e., the energy scale of inflation relative to the 
scale of new physics. As also noted in some previous works on the topic (e.g., [22]), the 
sensitivity to £ of a given probe is determined by two experimental parameters. One is 
the precision with which the power spectrum can be measured over the relevant range, 
specified by the error bars of the data and/or the width of the window functions. The 
other important parameter is the "lever arm", i.e., the span in fc-space over which the 
power spectrum can be measured, which is crucial for disentangling ripples in the power 
spectrum caused by trans-Planckian physics from a spectral tilt or an amplitude shift. 

Qualitatively, the precision of the data is the limiting factor in the case of very 
small £, while the lever arm determines, for a given e, the largest £ that can be detected 
uniquely as ripples in the primordial power spectrum. The limit imposed by the lever 
arm also means that the detectability of trans-Planckian ripples depends strongly on 
the magnitude of the first slow-roll parameter e. If e is too small, no detection window 
exists for £, as can be seen from equation (3.1), as well as more quantitatively from 
tables 1 to 3. In the same vein, if e = is consistent with data, then no bounds can be 
set on £. 

We have estimated the detection threshold of these ripples for data from future 
CMB experiments (Planck and a hypothetical cosmic variance limited experiment CVL), 
a galaxy redshift survey exemplified by the LSST, and finally a large scale 21 cm survey. 
Due to their very low noise and large volumes, both galaxy and 21 cm surveys will 
eventually reach a higher precision than CMB observations alone. This will allow us 
to reach down to smaller values of £ than ever before. However, large scale structure 
surveys will not be able to compete with CMB observations in terms of wavelength 
coverage. Therefore, in order to maximise the lever arm for the detection of trans- 
Planckian ripples, these surveys should be combined with CMB measurements. 

In an optimistic scenario with e = 0.01, we find that on the intermediate timescale, 
data from Planck+LSST will be open a detection window spanning 10~ 3 < £ < 0.02. 
In the more distant future the lower limit of the window can possibly be pushed down 
to 10~ 4 using data from a cosmic variance limited CMB experiment combined with 
the FFTT 21 cm survey. This same combination of data will also create a window of 
detection for e = 0.001 and possibly lower for the first time (10~ 4 < £ < 0.0025). 

So what are the implications for the scale of new physics A, which is presumably 
close to or at the Planck scale? Given that present data already constrains the scale of 
inflation to H < 2 x 10~ 5 Mp [44], we conclude that in the Danielsson model, there are 
no realistic prospects of ever detecting the traces of new physics in the primordial spectra 
if A > 0.2 Mp. For models in which the amplitude of the ripples is quadratic instead of 
linear in £, as predicted by purely effective field theory arguments in reference [18], the 
new physics would need to set in at energies smaller than ~2x 10 -3 Mpp 
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Nonetheless, we have shown that the fantastic precision of future CMB surveys 
like Planck, and large scale structure surveys such as LSST and FFTT will extend 
the detection window of trans-Planckian ripples by more than two orders of magnitude 
compared with present data. After all, the prospect of probing physics close to the 
Planck scale is a most tantalising one; a detection of trans-Planckian signatures would 
truly be one of the most important accomplishments of modern precision cosmology. 
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Appendix A. Details of the parameter error forecast 

The cosmological model used in our forecast consists of nine free parameters, the fiducial 
values of which are as follows: the baryon density Vt^h 2 = 0.02273, dark matter density 
Q c h 2 = 0.1099, Hubble parameter h = 0.72, redshift to reionisation z ve = 12, scalar 
spectral index n s = 0.96 and normalisation ln(10 10 y4 s ) = 3.18. The fiducial values of the 
slow-roll parameter e and the ripple parameter £ are given in the main text of the paper, 
while the fiducial phase <fi is 0. In addition, we use the following consistency relations 
on the tensor spectral index n t and normalisation A t : n t = — 2e, and A t /A s = 16e. 

We generate mock CMB data using the method discussed in reference [45]. The 
likelihood function £ is defined as 

*-max _____ I /"* I 

(A.l) 



r n-(Ci' 1 Ci)+ In JSl-n 

\Cp\ 



Xc 2 ff = -21n£=]T(2£+l)/ sky 

1=2 

where Cg = C^ v , (a,v = T, E, denotes the mock data covariance matrix, Ce = C% v ' + N^ V 
is the total covariance matrix comprising theoretical predictions of the CMB anisotropy 
spectrum Cg V and the noise power spectrum N^ u . The quantity n counts the number 
of observable modes. For observations in temperature and E-type polarisation, n = 2. 
Parameter estimation is performed by ^-minimisation, using a simulated annealing 
routine coupled to the Boltzmann code CAMB [46]. We often plot Ax^ ff (x) as functions 
of the parameters we wish to constrain (e.g., x — £, e). This is defined as 



A^ ff x) = -2ln 



£, 



(A.2) 



where £ max is the global maximum of the likelihood function £ defined in (A.l), and 
C^'(x) oc max C(x, y±, . . . , y^) (A. 3) 

is a projection of C(x, y±, . . . , y^) onto the one- dimensional subspace x by maximising 
along the yi, . . . ,yN directions. 

We define our "la", "2a" and "3a" intervals as the regions of x satisfying 
respectively A^cff < 1, 4, and 9. We emphasise that these intervals have no formal 
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probabilistic interpretation. However, if £W(x) is a Gaussian distribution and we 
assume flat priors on the model parameters, then these intervals are exactly identical 
to the 68%, 95%, and 99% credible intervals derived from a Bayesian analysis. See 
reference [47] for more details on Bayesian versus non-Bayesian interval construction. 
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